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Abstract 



We develop a theory of static BPS domain walls in stringy landscape and present a 
large family of BPS walls interpolating between different supersymmetric vacua. Examples 
include KKLT models, STU models, type IIB multiple flux vacua, and models with several 
Minkowski and AdS vacua. After the uplifting, some of the vacua become dS, whereas 
some others remain AdS. The near-BPS walls separating these vacua may be seen as 
bubble walls in the theory of vacuum decay. As an outcome of our investigation of the 
BPS walls, we found that the decay rate of dS vacua to a collapsing space with a negative 
vacuum energy can be quite large. The parts of space that experience a decay to a 
collapsing space, or to a Minkowski vacuum, never return back to dS space. The channels 
of irreversible vacuum decay serve as sinks for the probability flow. The existence of such 
sinks is a distinguishing feature of the landscape. We show that it strongly affects the 
probability distributions in string cosmology. 
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1 Introduction 



Soon after the invention of inflationary cosmology it was realized that inflation may divide our 
universe into many exponentially large domains corresponding to different metastable vacuum 
states pp. The total number of such vacuum states in string theory can be enormously large 
[21 El Ej • A combination of these two facts with the KKLT mechanism of vacuum stabilization 
|H] gave rise to what is now called the string landscape scenario [HJ. 

The purpose of this paper is to clarify some features of the landscape. In particular, one 
may wonder what are the properties of the domain walls separating different vacua, how these 
domains may form, and how large will be the fraction of the volume of the universe inside 
different domains. 

To this end, we will construct a broad class of static BPS domain wall solutions interpolating 
between different AdS or Minkowski vacua with unbroken supersymmetry. Quite generally, 
these BPS domain walls may interpolate between various minima, maxima, or saddle points of 
the scalar field potential. 

BPS domain wall solutions in 4-dimensional N = 1 supergravity were first studied in [7j. The 
methods used for that analysis have been later refined and generalized to a variety of models, 
for different number of dimensions and supersymmetries, see e.g. [SUSIE]]- Our present aim is 
to further extend these methods to make them suitable for investigation of the string landscape 
scenario. Following the logic of the recent developments, we will begin with an investigation 
of the supersymmetric vacua arising from flux compactifications prior to the KKLT uplifting, 
and later introduce the more physically relevant metastable de Sitter vacua. 

We start in Section 2 with an overview of some known features of domain walls in effec- 
tive d=4, N=l supergravity. The new results here include the derivation of BPS domain wall 
gradient flow equations in generic models with F-term potentials with arbitrary Kahler poten- 
tial, superpotential and number of moduli. We rewrite the gravitational action, including the 
boundary i^-terms in a BPS form: when the flow equations are satisfied, the integrand of the 
action defines the BPS domain wall tension. Afterward, in the spirit of the analogous work in 
five dimensions in the context of the AdS / CFT correspondence, where domain walls are viewed 
as dual to renormalization group flows [IT] , we present the BPS flow equations in the form of 
RG flows and provide the relevant "c-theorem" . 

In Section 3 we give examples of BPS domain walls where the potential has more than one 
critical point but there is no barrier between them. Our examples include an AdS minimum and 
an asymptotically Minkowski vacuum. Other examples have critical points at finite distance 
in the moduli space: an AdS saddle point and an AdS minimum, or an AdS maximum and 
an AdS saddle point. A particularly interesting case of flux vacua area codes in a Calabi-Yau 
hypersurface required the study of the flow of four moduli. They flow from one critical point of 
the potential, the Landau-Ginsburg fixed point, towards the second critical point with larger 
absolute value of the cosmological constant near the conifold. 
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In Section 4 we study a class of models in stringy landscape with a barrier between critical 
points. We have examples of a local Minkowski or AdS minimum and another absolute AdS 
minimum of larger absolute value of the cosmological constant. In case of AdS-AdS wall, 
we display the domain walls with critical and super-critical tension. The last case has the 
superpotential crossing zero between two AdS vacua. 

The results of our investigation in Sections 2-4 directly apply to the transitions between 
super symmetric AdS and Minkowski vacua. However, the most interesting parts of the string 
theory landscape correspond to dS vacua. These vacua can be obtained from supersymmetric 
AdS or Minkowski vacua by uplifting, which can be achieved, e.g., by adding some anti-D3 
branes |5]. This procedure breaks supersymmetry and may destabilize some of the vacua, in 
particular those which corresponded to maxima and saddle points. As a result, some of the BPS 
wall solutions may disappear after the uplifting, whereas some others will remain. However, 
instead of describing the static BPS domain walls, they will describe domain walls moving with 
acceleration. Such domain walls can be formed during the tunneling and bubble formation in 
the process of decay of metastable vacua, as discussed by Coleman and de Luccia (CDL) in 
[12]. We will discuss this process in Section 5. 

If the uplifting is relatively small, which may happen when we uplift a supersymmetric 
Minkowski vacuum, or a shallow AdS vacuum, as in the version of the KKLT scenario developed 
in Ref. then the properties of the domain walls do not differ much form the properties 
of the BPS walls prior to the uplifting. In this case we will be talking about near-BPS, or 
near-extremal walls and bubbles. The concept of near-BPS walls and bubbles is closely related 
to the concept of near-extremal black holes. Whereas at the extremal limit the BPS black 
holes are absolutely stable and have zero temperature, the near-extremal black holes have 
a non-vanishing temperature, they can evaporate and eventually approach a stable extremal 
limit. The physics of near-extremal black holes has been studied intensively over the years: a 
remarkable agreement with the physics of D-branes was established under condition of a small 
deviation from extremality |14j . 

In our new endeavor we will try to understand the rate of tunneling into a collapsing vacuum 
with a negative cosmological constant from uplifted slightly non-BPS vacua using the proximity 
of the large size bubbles to the BPS domain walls whose tension we can evaluate exactly. In the 
BPS limit the corresponding bubbles have an infinite radius, but the near-BPS ones have a large 
but finite radius [7j. The description of these bubbles can be performed using the deviation from 
the BPS limit as an expansion parameter. This simplifies investigation of bubble formation and 
vacuum decay in such models. In Section 5 we will find, in particular, that the vacuum decay 
rate of the uplifted supersymmetric Minkowski vacuum does not depend on the details of the 
model and is given by the universal expression e -s / 2 , where S = is the entropy of the 
dS state. This decay rate is much higher than the typical rate of decay of a dS vacuum to a 
Minkowski vacuum, which in many cases does not differ much from e~ s j^j- In application to 
our vacuum with Vds ~ 10~ 120 , we are talking about the difference in the decay rate by a factor 
~exp(10 120 /2). 
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The situation is especially interesting in a more general case when we uplift the models with 
several different AdS minima of the moduli potential. The AdS vacuum with a smaller value of 
the volume modulus experiences a bigger uplifting. The scale of supersymmetry breaking in the 
uplifted vacuum is related to the magnitude of the uplifting. As a result, the decay of the uplifted 
AdS vacuum to the collapsing space with a negative vacuum energy usually occurs much faster 
than the typical dS rate of the uplifted Minkowski vacuum ~ e~ s / 2 . Under some conditions, 
the decay rate can be quite large, with the suppression factor e~ s = exp(— ^j) ~ e - 10l2 ° 

being replaced by exp 

(~rfeS)- Here C = °W and \ VAi3 \ is the depth of the AdS 

minimum 

prior to the uplifting, which may be a hundred orders of magnitude greater than Vds ~ 1CT 120 . 
In the KKLT-type models, the depth of the AdS vacuum prior to the uplifting is of the same 
order as fny 2 , where 7723/2 is the gravitino mass after the uplifting [13] • Therefore in the 
class of models studied in Section 5, the vacuum decay rate due to the tunneling to a collapsing 
universe is related to the gravitino mass, P ~ exp(— 0(m~^ 2 )^ . After the tunneling, the universe 
within each bubble with negative vacuum energy collapses within a microscopically small time 
~ I VAds]^ 1 ^ 2 - The condition that the lifetime of the uplifted vacuum is greater than the age of 
our part of the universe may serve as a vacuum superselection rule. 

An important feature of the decay of dS space to a Minkowski vacuum or to a collapsing 
universe with a negative vacuum energy density is that the part of space that experiences such a 
decay never returns back to dS space. This is quite different from the back-and-forth transitions 
between various dS vacua, which may keep the system in a kind of thermal equilibrium. In this 
respect, the channels of irreversible decay of an uplifted dS vacuum to a Minkowski vacuum or 
to a collapsing universe serve as sinks for the probability flow in the landscape. In Section 6 
we will show that the existence of such sinks has important implications for the calculations of 
the probabilities in the landscape. 

In particular, if one ignores the different rate of expansion of the universe in different vacua, 
i.e. uses a comoving volume distribution, then in the absence of the sinks, the fraction of 
the comoving volume of the universe in the state with a given vacuum energy density V > 
asymptotically approaches e ~ e - ^. This distribution has a sharp maximum at the smallest 
possible V. Meanwhile if the vacuum with the smallest possible V is not the ground state but 
a sink, and the probability of tunneling to the sink is sufficiently large, then the fraction of the 
volume in the vacuum with the smallest V will be exponentially suppressed, and the maximum 
of the probability distribution will be at larger values of V. 

Walls and bubbles considered in our paper interpolate between domains with different values 
of the scalar fields, but with the same values of fluxes. In Section 7 we will briefly review our 
results and discuss a possibility of their generalization for the transitions between vacua with 
different fluxes. 
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2 Supersymmetric Domain Walls in N = 1 supergravity 



2.1 Preliminaries 

In this paper we are going to analyze the domain wall solutions interpolating between different 
supersymmetric vacua of the landscape of string theory flux compactifications to 4 dimensions. 
We will focus on the case of N — 1 supergravity. BPS domain wall solutions with 1/2 of 
unbroken N = 1 supersymmetry in the minimal 4-dimensional supergravity have been discussed 
in [7j. The equations presented there refer to the case of a single scalar field. Here we are going 
to extend this analysis to a more general setup allowing for the possibility of many moduli fields 
with generic metric in the moduli space. Of course, in any sufficiently small patch, there is a 
possible reparametrization of the moduli space so that there is a single scalar field supporting 
the domain wall. In a rather generic setup, however, this is not possible globally, and therefore 
the equations we will propose in the following are more useful in this respect. The analysis 
presented here is mainly based on the approach presented in |S] for the domain wall solutions 
of the minimal 5-dimensional supergravity. 

We briefly recall the main ingredients needed in order to discuss supersymmetric bosonic 
configurations of N=l supergravity coupled to chiral multiplets. The bosonic sector of the 
action is 

S = J d A x^~g Qi? - g^z'd^z? - V{z, z)J , (2.1) 

where the potential can be expressed in terms of a superpotential W(z) as 

V = e K (g iJ DiWD]W -3\W\ 2 ), (2.2) 

or, in terms of the covariantly holomorphic superpotential Z = e K ' 2 W, as 

V = \DiZ\ 2 -3\Z\ 2 . (2.3) 

The covariant derivative is defined as Di = <9, + diK on holomorphic quantities (like W) and 
as Di = di + \diK on covariantly holomorphic quantities (like Z). We will not consider vector 
multiplets and possible D-term extensions here. 

The relevant supersymmetry transformations (of the gravitino and modulini) are 

5%l = ^ + \ujf lab + l -A^ L + ^ 2 W^e R , (2.4) 

% = \fd^e R - ^e K / 2 D]We L , (2.5) 
where = | (d^djK — d^z l diK) is the £7(1) Kahler connection. 
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2.2 The BPS equations 

In the following we will consider flat 1 space-like domain wall solutions of the supersymmetry 
conditions and of the supergravity equations of motion. For this reason we are going to use the 
following metric Ansatz 

ds 2 = a 2 (r)(-dt 2 + dx 2 + dy 2 ) + dr 2 . (2.6) 

Given the domain wall structure, we split the indices \x = {i,r}, a = {i, r}. For the metric 
(|2.6jl there is only one non-trivial component of the spin connection 

uf = a!5], (2.7) 

where leading also to the Ricci scalar (which will become useful later) 



R 



,2 // 

a \ a 



(2i 



Plugging ()2.6|) into the supersymmetry equations ()2.4j) - (j2.5|) and assuming that also the scalars 
z % depend only on the radial coordinate r, one obtains 

HrL = d r e L + l -Afe L + l -e K ' 2 W lT e R (2.9) 
6i/>iL = ^IHrjL + le^W^eR (2.10) 

5xl = ±Y^(r)e R -±e K ' 2 gWjWe L , (2.11) 

where di€i = was assumed. Using also 7, = -7$, one can start solving the various equations. 
From (|2.1(J|) . splitting the complex superpotential into norm and phase 

W = \W\e ie{r) , (2.12) 

we get the projector condition 

lrJ L = Te ie e R . (2.13) 

from which the BPS equation follows 



a' 



- = ±\W\e R / 2 = ±\Z\. (2.14) 
a 

The condition ()2.13|) can be assembled into a real projector V (satisfying V 2 = V) acting on 
the full spinor as Ve = e: 

V = ±— ^-r \lr ± (COS 6 + Z 75 sin 6)] . (2.15) 

2 cos ~ 



1 Curved domain wall solutions have been considered for instance in PITS], 
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The radial equation on the gravitino gives a differential equation on the spinor parameter 

(2.16) 



-\A*e L ±\e K l*\W\e L . 



Here, the projector condition was used. Moreover, one can use the projector condition in- 
side ([2.160 . or take the charge conjugate (using the Majorana condition). The two resulting 
equations are 



-R 



-i9'e R - l -A?e R ±h K / 2 \W\e R , 
~A*e R ±\e K '*\W\e Rl 



whose consistency implies 



& = -A 



B 



(2.17) 
(2.18) 

(2.19) 



The solution is given in terms of an arbitrary constant left-handed spinor parameter eo 
(hence the 1/2 supersymmetry preserved): 



(2.20) 



The BPS equation (|2.14j) has been used to express the integration over the superpotential in 
terms of the warp factor and ([2.190 for the integration over . In the last formula we have 
used the definition of the charged conjugate spinor e$ that follows from choosing the charge 
conjugation matrix C = jrfl, with (3 = 270. In this way the spinor e satisfies the relation ([2.130 . 

Finally, using again the projector condition in the modulini equations, one finds the BPS 
equation for the scalar fields 

z i> = Te K/i+M g iTD~w. (2.21) 

We can now also remove the dependence on the phase 9, by using the holomorphicity of the 
superpotential: 

dW = B\W\ + i\W\B9 = 0. (2.22) 

This implies that the BPS equations for the scalars can take the form: 

z i> = Te K ^ g ^{2dj\W\ + djK\W\) = T2g% (e K/2 \W\) = ^2g%\Z\. (2.23) 

Summarizing, the supersymmetry conditions that have to be satisfied in order to specify the 
domain wall solutions are 



(2.24) 



The other consistency condition we derived from the solution of the supersymmetry equations, 
namely ([2.190 . is identically satisfied. This can be seen by using the fact that 9 is the phase 




7 



of the superpotential, and by equating with (|2.19jl the condition coming from its definition, 
9' — ~ {^—^y + ^pfj- The resulting equation is an identity upon using (|2.24j) . 

Domain walls m. N — 1 supergravity were studied before in [7] in a more restrictive setting, 
so we compare our results. First there is a change in the signature. This implies that one should 
change the sign for every upper index and put an i for each gamma matrix in the supersymmetry 
equations. This can be re-absorbed in sending Wcv —inhere- Then there is a difference in 
the metric. On the solution they have A = B = a 2 , but also their radial coordinate z is related 
to ours by ^ = vA Also, the decomposition of the superpotential in norm and phase implies 
that 6cv = Altogether the equations in [7j become 

& KI2 \W\, (2.25) 

-(e K / 2+w g TT D^W, (2.26) 

l - (r'd T K - T'dfKj . (2.27) 

It is easy to see that these equations are equivalent to the ones presented above, when restricted 
to a single scalar field flowing. 



of_ 
a 

r = 

6' = 



2.3 A "c-theorem" for N=l/2 domain walls 

We find it useful to present the domain wall equations in the form of renormalization group 
equations. For this purpose we will describe the flow not as a function of r coordinate but as 
a function of the warp factor. We multiply Eq. (j2.21|) on 4 and find, using the Eq. (j2.14|) 



V = -g«S^. (2.28) 
a' W 



and 

dz 



a— = —g v dj(K + In W) 



(2.29) 



We will find below that in certain cases the domain wall BPS equations in the form ()2.24j) are 
useful. Sometimes we will use equations in the form ()2.29j) . particularly in the case of many 
moduli participating in the flow and W not crossing zero. 

A careful inspection of the BPS conditions ()2.24|) shows that we can associate to the domain 
wall solutions a monotonic function, which therefore specifies the direction of the flow. As it 
is known since the analysis of the 5-dimensional domain walls and the interpretation of such 
solutions as renormalization group flows in the context of the AdS / CFT correspondence, there 
is a c-function that is given by the first derivative of the logarithm of the warp factor c = (log a) 1 
[TT] . By using the BPS equations we easily see that 

d = = ±\Z\ = ±{z i 'd i + fdf)\Z\ = Ag i ^d l \Z\dj\Z\ > 0, (2.30) 
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where the last inequality follows from the properties of the scalar metric. 

There is, however, another interesting quantity that in most cases works as a c-function. 
This is the square of the covariantly holomorphic superpotential, which is related to the square 
of the mass of the gravitino at a generic point of the moduli space 

\Z\ 2 = e K \W\ 2 = ±M 2 3/2 (z,z). (2.31) 

We can now look at the variations of \Z\ 2 along the flow. In order to compute it, we actually 
evaluate (|^| 2 )' = z 1 ' di\Z\ 2 + z 1 ' dj\Z\ 2 and use the BPS equation (|2.23p . After some easy algebra 
we get that 

(\Z\ 2 )' = ^8n\Z\g is di\Z\d^Z\. (2.32) 

Due to the positive defmiteness of the metric of the scalar fields, this expression is clearly 
either always increasing, or always decreasing along the flow, depending on the sign of the BPS 
equation, unless the flow crosses points where the superpotential vanishes. For these special 
cases the BPS equations change sign as we cross the zeros of the superpotential and therefore 
\Z\ 2 is first decreasing until it reaches zero and then starts increasing again. 



2.4 BPS form of the action 



The supersymmetric domain walls studied above are static 1/2 BPS configurations of N — 1 
supergravity. In an analogous situation for static 1/2 BPS cosmic strings in N=l supergravity 
[To] it was possible to rewrite the gravitational action, including the boundary terms, in the BPS 
form. The action therefore was given by some squares of the BPS conditions and total derivative 
terms. A saturation of the BPS bound is equivalent to solving first order differential equations 
for unbroken supersymmetry. The on shell action therefore is equal, up to a sign, to the total 
energy of the string. Therefore the tension of the string, the energy per unit area, can easily be 
derived from the BPS action. The metric of the cosmic string is ds 2 = —dt 2 +dz 2 +dr 2 +C 2 (r)d8 2 
and the relation between the —S on shell and a tension \x is 

-5~ trinSS = Jdtdz J drd6^T° = J dtdzfi, (2.33) 

where \l = 27r£ and £ is a FI term in the action. 

For flat domain walls the metric is ds 2 = a 2 (r)(—dt 2 + dx 2 + dy 2 ) +dr 2 . Therefore one would 
expect the following relation between the action and the tension (the energy per unit of area) 

-SoSi waUs = J dtdxdya 3 (r)a = Jdtdidya. (2.34) 

In the case of cosmic strings t and z have flat metric, therefore the volume of the area is 
just J dtdx. However, for our domain walls we have to distinguish the physical volume 
from coordinate volume and we have to ignore (divide by) the physical volume which is 
J dtdxdya 3 (r) = J dtdxdy. 
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In [Z] the tension of the BPS domain walls was calculated using a rather complicated Nester 
construction. The formula for the tension/energy density of the wall between two AdS critical 
points is given for two cases: in one case the superpotential does not vanish between critical 
points, in the other case it vanishes at some intermediate point: 



2 AdS vacua, W ^ in between 



V5v 



a : ^,|A|r-|A| 



1/2 



(2.35) 



2 AdS vacua, W = in between 



a = 4flA|i /2 + IA 



1/2 



x/3 V 



;2.36) 



Here we show that for domain walls we can put the supergravity action in the standard BPS 
form, namely as the integral of the square of the BPS conditions plus boundary terms. We 
start with the gravitational action supplemented by the Gaussian curvature at the boundaries. 



S 



cf x l-R - grp^d^z 3 - V(z, z) 
+ J d 3 x (y/— det h K\ r=+cx> - V— det h K\ 



(2.37) 



Here h is the metric at the boundaries. In our case the boundaries are at r = const surface. 

yy = n u = a ~ 2 and K ab = ln^^ 



Therefore a/— det h = a 3 and K = g ab K a b where g xx = g 
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and = -^=y^j . One starts by plugging (|2.8|) in the Einstein term and using that all the 
quantities depend only on the radial coordinate. Then, we try to re-combine the various terms 
in the form of the BPS equations (|2.14|) - (j2.21|) plus some boundary term. The action (J2.37|) 
becomes 



S 



d 3 xdr 



-a 3 I 3 ( -V + 3- I - a 3 giJ z l/ z 3 ' + 3a 3 \e K / 2 W\ 2 - a 3 e K g i W l WD~W 



a 



+ I d 3 x(a 3 )'\ r=+00 - (a 3 ) 



3VI 



(2.38) 



This can now be recast in the form of BPS equations squared plus additional boundary terms 



as 



S 



d xdr 



3a d 



± e R / 2 \W\ - a 3 9lJ U' T e K / 2+ie g ik DW) T e K/2 ~ lt 'g* D k W) 



d 3 xdr [2(a 3 )'\e K/2 W\ + a 3 z'' d^^W)^ + a 3 z%(e K/2 W)e ie + a 3 {e K/2 )'\W\] (2.39) 



10 



Finally, the second line becomes a total derivative by using fl2.22j) . and the total action reads 



S 



d xdr 



3cT ( — ± \Z\ 

a 



gi Jz l, T2g fk d- k \Z\) (z^2g> k d k \Z\) 



T J d 3 xdr^-[a 3 (2\Z 



[2.4Q) 



where we recovered the BPS equations in the first line and we have a boundary term contribution 
in the second line. One can define the energy density functional a, since we have a static 
solution, so that S\ ons h e ii = — f d 3 xa 3 a. If the superpotential does not vanish between the 
critical points we find 



d x a a 



d 3 x dr 



-3a 3 \ — + \Z\ 



z Ll -2g ik d k \Z 



+ / d 3 x{a 3 2\Z\ 



r=+oo @> ^\^\r=—oo 



)■ 



(2-41) 



In case that \Z\, 



r=+oo 



-\z\ 



> we have to make a choice of the sign in BPS equations such 



that the tension is positive, as shown in Eq. (j2.41|) . Thus when \Z\ does not vanish between 
fixed points and the BPS equation of motion are satisfied with the upper sign the energy density 
of the wall is given by 



a 



2{\Z\ 



\Z\ 



(2.42) 



Note that A = 3\Z\ 2 . So our answer is exactly as in eqs. ()2.35|) . in agreement with |7 • It was 
explained in [7j that the BPS domain walls with this tension are related to CDL bubbles with 
infinite radius. 

Now we assume that there is a solution of the BPS equations where the superpotential 
vanishes at some point ro between the AdS critical points. The boundary integral in Eq. ()2.40j) 
must be split into two terms since one part of the flow proceeds with one sign in the BPS 
equations and for the other part, after Z = point, the sign is flipped. Therefore the answer 
for the energy density in this case is 



a 



Z(r )=0 



(2.43) 



Noting that A = 3|2| 2 , we have derived the correct tension as in Eq. (J2.36)) in agreement with 
[7j. This case, according to 7j is not related to CDL tunneling. 

We would like to conclude this subsection by the observation that our new derivation of the 
BPS domain wall tension follows the pattern discovered for the BPS cosmic string in ^H]- We 
find the derivation using the BPS form of the action significantly more elegant than using the 
Nester construction. 
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3 BPS domain walls between vacua not separated by 
barriers 



In this section we are going to give some examples of supersymmetric domain wall solutions 
for some of the generic models which capture the most important features of the landscape 
of flux vacua. We start with the simple scenario with just one Kahler modulus and all the 
complex structure moduli frozen: the KKLT model j^]. This model has one supersymmetric 
AdS vacuum at a finite point of the volume modulus and a Minkowski vacuum at infinite 
value of the volume modulus. We will solve the BPS equations numerically and display the 
corresponding domain wall. After this simple example we show a more involved setup where also 
the complex structure moduli and the axion-dilaton are allowed to flow: the STU models. These 
include a wide variety of models (and also the KKLT model) as a special case. The case which 
we will solve numerically has two different AdS critical points at finite distance in the moduli 
space and three scalars are flowing between the critical points, two of them proportional to each 
other. Finally we focus on a model where the effect of the Kahler deformations is neglected, 
but the form of the complex-structure deformations is more involved. The flux vacua in such 
models are described by a simple set of "new attractor equations" [TZ| relating the points in 
the moduli space where the moduli are stabilized by fluxes 2 . These flux vacua attractors have 
some resemblance as well as some differences with the black hole attractor equations [20J. We 
will focus here on a model of a Calabi-Yau threefold obtained as a hypersurface in WP^ X 112 . 
This model has two different AdS critical points [21 at finite distance in the moduli space. We 
will find out how four scalars in this model are flowing from the Landau-Ginsburg AdS critical 
point towards the near conifold region of the moduli space. 



3.1 The simplest KKLT model 

The KKLT model jS] suggests a way to generate meta-stable de Sitter vacua in string theory 
starting from supersymmetric Anti de Sitter minima with all moduli stabilized due to fluxes 
and using non-perturbative effects, such as gaugino condensation and/or D3-instantons. The 
uplifting of these stabilized AdS vacua to dS vacua is possible under certain conditions by 
the addition of positive energy density sources, like anti-D3 branes. The simplest realization 
has a single volume modulus p = e^^ + ia. After all the complex structure moduli have 
been fixed by the fluxes, the superpotential for the volume modulus (whose Kahler potential is 
K — — 3 log(p + p)) is 

W = W + Ae~ ap , (3.1) 
with Wo, A and a parameters depending on the details of the setup. 

The corresponding potential, before uplifting, has one AdS supersymmetric critical point 
at finite volume and an asymptotic Minkowski vacuum at infinity of the volume. It is easy to 

2 These new attractor equations for IIB vacua have also been extended to the case of non-Kahler compacti- 
fications of the heterotic theory |18|. Similar extensions for the black hole case have been discussed in |19|. 
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check that we can consistently set the axion a to zero, since this is a fixed point for any value 
of the volume <fi. For the sake of simplicity we will also take the parameters A, a and Wo to 
be real (and positive the first two A = 1, a = 0.1, negative the last one Wq = —0.0001). The 
potential for this case and the covariantly holomorphic superpotential Z are given by Fig. 13.11 
as functions of the volume modulus. 



V 




Figure 1: KKLT potential V (blue upper curve) and Z (red curve below) as functions of the modulus. 
One can see that the potential has an AdS minimum at <j> = <fi* ~ 4.09854 and tends to Minkowski 
limit at large <p. Z is always negative in the interval between AdS and Minkowski vacua. Potential 
energy density V is shown in units 10 -15 Mp, whereas the covariantly holomorphic superpotential Z 
is shown in units 10 _8 Mp. 



One finds that for these parameters during the total flow the covariantly holomorphic su- 
perpotential Z is negative and approaches zero at infinity. We can reduce the problem to a 
one-dimensional gradient-flow, where the BPS equations follow by using 

Z\ ^ = -\Z\ = —=- (w + Ae- aeV * . (3.2) 



We solve numerically the flow equations in the opposite direction (for technical reasons), starting 
with the AdS minimum towards the infinitely far Minkowski. For this purpose we have to use 
(J221: 

a\r) = -a(r)\Z(r)\, (3.3) 
0'(r) = d^\Z(r)\. (3.4) 

The AdS critical point is at <fi = <ft* ~ 4.09854, and for values slightly above <p* we get a domain 
wall solution that drives — ► oo and Z — > 0. In this case we obviously have just one basin 
of attraction and the attractor point is the supersymmetric AdS vacuum. Asymptotically, at 
large r, the value of <fi grows, and \Z\ becomes 

|Z|~ce-^, c=\W \/V8>0. (3.5) 
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We can therefore solve the BPS equations above in this limit by 



a ~ a {r -r )3. 



-^log(3c(r -r )). 



It should be noted that 



1 



(3.6) 



(3.7) 



a 3(r — r ) ' 

which means that the curvature of the space-time (|2.8j) vanishes asymptotically, as it should be 
in Minkowski space. 




r 




R 



r 



10 12 14 



r 



10 12 14 



Figure 2: The KKLT wall interpolating between the AdS and flat Minkowski space which are at 
infinite distance from each other in the moduli space. We plot the scalar <j>, the covariantly holomorphic 
superpotential Z (in units 10~ 9 ), the warp factor a and the curvature R (in units 10~ 15 ), all as functions 
of r. Here r is the coordinate of the domain wall configuration given in Eq. (|2.6|) : we show r in units 
10 7 . 



In Fig. El we plot the wall: first on the top right we have a scalar interpolating between 
the AdS value = 0*^ 4.09854 and infinity. On the top right we plot Z as the function 
of wall coordinate r. It is always negative but here we see clearly that it takes a constant 
negative value at AdS and starts interpolating towards zero at Minkowski vacuum. The next 
plot on the bottom left shows the warp factor a(r). Finally on the right bottom we show the 
curvature of the domain wall: it is vanishing at large values of modulus in the flat space and 
grows towards the decreasing values of modulus and sharply drops down to its negative AdS 
value upon reaching the AdS vacuum. 

From this example we learn that in the absence of a barrier between the AdS and Minkowski 
vacua there exists a static domain wall where the flow of the volume modulus from the infinite 
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value at Minkowski towards the finite value at AdS is controlled by 1/2 of unbroken supersym- 
metry. The tension of this wall is equal to a = -^Al 1 / 2 where A is the cosmological constant 
of the AdS minimum. 



3.2 STU models 

The previous example had only one AdS vacuum. However, in the landscape of flux compactifi- 
cations there are also examples with multiple AdS vacua and therefore it is interesting to display 
an example with multiple basins of attraction to different AdS vacua and the supersymmetric 
domain wall solutions interpolating between them. In order to find a workable example, we 
need an analytic superpotential admitting multiple AdS vacua. Such an example is provided 
by various models of toroidal orbifolds analyzed in [22] • The generic moduli content is given 
by the axio/dilaton S and by the complex-structure moduli U l and 3 Kahler moduli T l , with 
Kahler potential 

h (2,l) 

3 untwisted 

K = -\og(S + S)-Y,log(T i + T)- ^(W + TT). (3.8) 

i=l i=l 

For the sake of simplicity we focus on the case of one untwisted complex structure modulus, as 
for instance in the Z2 x Z3, Z2 x 7*$ orbifolds or in the 1iq-u orbifold with SU(2) 2 x SU(3) x G2, 
SU(3) x SO(8), SU(2) 2 x SU(3) x SU(3\ or 517(2) x SU(6) lattice. The superpotential is a 
polynomial in the complex-structure and dilaton moduli that comes from the fluxes, added by 
some non-perturbative exponential terms involving the Kahler-structure moduli. If we consider 
the cases with just one complex structure modulus, the full untwisted superpotential reads [22] 

3 

W = a + bU + cS + dSU + ^2g i e- hiTi . (3.9) 

i=i 

Here the parameters a, b, c and d are related to the fluxes (they are the coefficients of the 
real sections of the 3-cycles), while the parameters and hi are to be associated with non- 
perturbative effects. This potential includes the one-modulus KKLT scenario as a special 
example. 

We can then look for supersymmetric vacua, satisfying DsW = 0, DyW = and D T iW = 0. 
For simplicity, let us consider the case that g\ = g, hi = h and look for vacua where T % = T. 
Then, the supersymmetry conditions become 

c + dU = j^W, (3.10) 
b + dS = jj^W, (3.11) 
-ghe~ hT = j^W. (3.12) 
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A generic solution for W 7^ of these equations requires that 



U=^S = kS, (3.13) 



and if we further look for the case where ImT = ImS = 0, then we also obtain that 



ReT = -- log 
h 



a-kdRe S 2 
35 



(3.14) 



3cReS(b + dReS) 
Ke " h(ab-cdReS 2 ) ' ld 5j 

When the values of the fluxes are chosen appropriately, one can obtain 2 distinct vacua at 
values of s = ReS > and t = ReT > 0. One instance is the choice a = 1, b = 2, c = 
4,d = -1,5 = -10, /i = 1/100. The critical points are at {s = 0.496855, t = 299.992} and 
{s = 0.749598, t = 264.799}. (The value of the size of the manifold is tuned by 1/h. The 
smaller h, the bigger t. Increasing g on the other hand allows for smaller values of s.) 

We can then use the BPS equations to obtain the interpolating solution between them. 
Having seen how the critical points are obtained, we would like to reduce the problem to a 
two-dimensional system depending only on s = ReS, t = ReT. Both critical points lie on the 
section of the moduli space where all the axions are set to zero. This condition can be kept 
consistently all along the flow. When the axions are vanishing, W = \W\ and therefore 8 = 0. 
This condition can be preserved along the flow as 9' = if the axions are vanishing. This also 
implies, since the metric is diagonal, that the real and imaginary parts decouple consistently. 
Given the symmetry of the potential, we can also assume that T % = T all along the flow. The 
BPS equations therefore reduce to 

s' = T8s 2 d s \Z\(s,t,u), 

t' = ^8t 2 d t \Z\(s,t,u), (3.16) 
u = ^f8u 2 d u \Z\(s, t, u), 

and, by using (|3.13jh one can prove that the differential equations for u and s are proportional 
to each other. This means that we can further reduce the problem to the description of the 
gradient flows for the flow potential where s = e^^ 1 and t = e^^ 2 . In these variables, 

the BPS equations become 

0', = Td^Zm 

a' ( 3 - 17 ) 

- = ±1*1(0*) 

a 

where the absolute value of the central charge reads 



\Z\(<f,i) = \ / J> >/8(-*i-3* 2 ) f a + £ e ^(26 + de^) + 3ge h ^ 2 ) . (3.18) 

V o^C \ / 

Notice the double exponentials due to the non-perturbative terms in the Kahler moduli. 
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Given this flow potential, we can look for the special case where we have two distinct critical 
points, for instance when the parameters are chosen as above: a = 1, b — 2, c = 4, d = —1, g = 
— 10, h = 1/100. For these parameter values (jH.18|) shows two supersymmetric AdS critical 
points with different values of the cosmological constant and an asymptotically Minkowski 
vacuum. 

Although the BPS equations can be written in a closed form, the integration can only be 
done numerically. In Fig. |3]we show the contour plot of Z, where the (f)\ modulus (proportional 
to the complex structure and to the dilaton s = e^^ 1 ) is on the horizontal axis and 02, 
proportional to the Kahler modulus t = e^^ 2 is on the vertical one. We plot the total potential 



2 .06 
2 . 04 

2 . 02 
2 

1 . 98 
1 . 96 



-0.4 -0.3 -0.2 -0.1 

Figure 3: Here is the contour plot of the covariantly holomorphic superpotential iT(0i,02) and the 
gradient flow between the two critical points. i^(0i,02) is everywhere positive between AdS critical 
points, it vanishes only at the asymptotically Minkowski vacuum. 

on Fig. |U We can clearly see from the pictures that there are two distinct critical points, and we 
have plotted the numerical solution interpolating between them. As it is clear from the picture, 
the two critical points have different behavior. Using the nomenclature coming from 5d, one is 
approached by the flow in an ultra-violet direction (i.e. the warp factor goes as a ~ e r —>■ +oo, 
for r — > oo), the other in an infra-red (i.e. the warp factor goes as a ~ e r ~ 1 for r — > — oo). An 
ultra-violet critical point is a minimum of the absolute value of the superpotential (minimum 
with positive W or maximum with negative W) , an infra-red critical point is a maximum. This 
implies that there is no need for the superpotential to cross zero. The warp factor indeed shows 
no critical points, as it always increases along the flow. The superpotential will necessarily cross 
zero along the flow and therefore the BPS equations will change sign only when the interpolation 
must be done between two infra-red (two ultra-violet critical points are excluded as the wall 
should be singular to interpolate between them). This is the only way to obtain a smooth 




17 



V 



-6.31 



-6.28 



-6.29 



20000 




80000 



100000 



r 



Figure 4: Flow of the potential V(<f>i(r), <fa,{r)). 



domain wall interpolating between them. A similar phenomenon occurred in the realization 
of the Randall-Sundrum scenario with a thick brane in (231- We will have examples of such 
domain walls in Sec. 14.2.21 

From Fig.Elwe can also clearly see the two basins of attraction for the AdS and Dine-Seiberg 
Minkowski critical points and the AdS repeller (with a marginal line of attraction separating 
the two basins). 

This shows that there is one main basin of attraction for the AdS vacuum with the larger 
absolute value of the cosmological constant. There is a line of attraction for the other critical 
point, which is really a saddle point of the flow potential. Finally there is another basin 
of attraction for the values of the moduli that leads to the decompactification of the space. 
The size modulus t runs to infinity at the same time as the dilaton goes to zero, so we get 
a decompactification of the space in the weak coupling regime. This example shows that the 
landscape attractors are different from the black-hole attractors. For the latter, supersymmetry 
vacua are always minima of the superpotential W and therefore one always has clear distinct 
basins of attractions and there are no repeller-like vacua. Moreover, if one is looking at the 
5-dimensional case there is always only one basin of attraction for the vector-multiplet moduli. 

3.3 Example of a Calabi— Yau hypersurface 

The panorama of the landscape can be further refined by considering a representative example 
in the class of type II compactifications on Calabi- Yau threefolds, defined as hypersurfaces in a 
weighted projective space. When fluxes are active for the RR and NS three-form field strengths 
on the orientifold of these hypersurfaces, a non-trivial superpotential is induced for the dilaton 
and complex structure moduli, leading to the many interesting physical properties that have 
driven a massive research activity (for a review and bibliography see 

These models are particularly relevant even for a single deformation parameter because of 
the results obtained in showing that, quite unexpectedly, the density of flux vacua is not 
uniform in moduli space, but rather it peaks around the conifold singularity, x/j — 1. Indeed, 
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upon extensive numerical computations, it has been shown that the conifold region plays the 
role of "attractor" basin and it serves as accumulation region for flux vacua. This supports the 
analytic predictions made in |2E] • On the other hand, also the Landau-Ginzburg singularity at 
-0 = enjoys very peculiar properties, deriving from it being a fixed point of the modular group. 
In general, it is the special point in moduli space where the number of vacua with vanishing 
superpotential W = and discrete symmetries can be of the same order as the total number of 
vacua (when non-zero) [21]. For our purposes, it will be convenient for computations because 
of our need to use perturbative expansions for the CY periods, that are known around this 
point [2E]- Thus, in spite of neglecting the Kahler structure moduli that in principle would be 
present in these systems, this case is particularly interesting as a representative of the large 
amount of possible flux vacua that are making the landscape so rich. 

Another important related feature is that one may encounter multiple attraction basins 
[2*Tj . These consist of distinct supersymmetric vacua that can be found within a finite region 
in moduli space, or "area code" , for fixed values of the F( 3 ) and if (3) fluxes. We thus proceed 
to identify such a situation and to construct the domain wall solution interpolating between 
different area codes. As a start, we quickly recall some necessary formulas, while referring for 
instance to |2Z| for the relevant lore and more conventions. 

One examines in general the ^(2,1) complex structure deformations. As usual, one chooses a 
symplectic basis {A a , B b } for the 63 = 2/i( 2 ,i) + 2 homology cycles (a, b = 1, ... , /i( 2 ,i) + 1), with 
dual cohomology elements {a a , /3 b } such that 



a b = 5 a b , p a = -5 a b , a a A(3 b = 5 b a , (3.19) 

A a J B b J 

and one assembles into a symplectic vector the periods of the unique holomorphic three-form 
ft, 



X a = / ft, T h = / ft, Tl(X) = (F a ,X a ). (3.20) 

J A a JB b 

The X a are projective coordinates on the complex structure moduli space and T b = dbJ-'(X), 
where F{X) is the holomorphic prepotential. 

The Kahler potential for the complex-structure and axion-dilaton moduli space is given by 
K = - log (i J ft A ft J - log(-i(r - f )) = - log {-iYl ] • £ • II) — log (-i(r - f)) , (3.21) 

where r is the axion-dilaton and X = ^ ^ y su P er P°tential induced by non-trivial 

3-form fluxes reads 

W= [ G (3 )(r)Aft(X) = (27r)V(/-r/ i )n(X), (3.22) 

J M 

where G( 3 ) = — tH^ and (/, h) are the integer valued flux components 

F (3) = -(2n) 2 a'(f a a a + f a+h2 , 1+ M , #(3) = -(2tt) V(/i a a a + h a+hxi+1 {3 a ) . (3.23) 
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We will take for convenience (2k) 2 a' = 1. 

Let us consider in the following M to be the Calabi-Yau locus defined by the equation 



} j x\ + 2%% - 6lf) X^X^X^X^ = 0, Xi e WP^ l l l 2 . 



(3.24) 



i=l 



The supersymmetric vacua are found by solving the "F-flatness" condition, that is a (Kahler 
covariantized) extremization with respect to r and the single complex-structure modulus ip, 

D T W = (d T W + d T KW) = , D a W = (d a W + d a KW) = , (3.25) 

leading to the equivalent conditions 

(/ - fh) ■ u(x) = (/ - rh) ■ (^n + a^n) = o. (3.26) 

Alternatively, one can solve the flux vacua attractor equations [Fj\. One finds that the moduli 
space of this model presents Landau-Ginzburg, conifold and large complex structure singular 
points respectively at if) = 0, ip — 1 an d if) — > oo j2Tj . 

The monodromy group, T, of the complex structure moduli space has two generators: A, 
which generates phase rotations if> aif> with a = exp(2ni/6) around the Landau-Ginzburg 
point at if) = 0, and T which corresponds to the logarithmic monodromy JF 2 — > JF 2 + X 2 around 
the conifold singularity if) = 1. By itself, A generates a Z 6 C T subgroup, with an associated 
fixed point at if) = 0; T, on the other hand, is of infinite order. The monodromy matrix A 
generates rotations by a root of unity around if) = 0: 

AU(ip) = aU(aif)) 

and is explicitly given by 

/ 1 -1 1 \ 
10-1 
-3-3 1 3 
-6 4 1 -3 ) 

The periods admit, for \i\)\ < 1, an expansion in a Picard-Fuchs basis 



.4 



V 



(3.27) 



(3.28) 



Wi(lf)) 



(2m) 



E 



exp(^pn)r( 



6 ^ r(n)r(i - |) 3 r(i - f ) 



6a l 
2V3 



if) 



n-l 



where a is the 6 th root of unity 



a 



exp 



2ni 

IT 



(3.29) 



(3.30) 



The usual choice for the periods in this basis appropriate for the polynomial M is w T = 
(w2(ip),wi(ip),wo(ip),w^(ip)). Then, the periods in the symplectic basis n(X) are related to 
the Picard-Fuchs periods by a linear transformation 



n 



?2 
X 1 

\X 2 J 



m ■ 



\w 5 J 



(3.31) 
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The matrix of transformation from the Picard-Fuchs to the symplectic basis has rational ele- 
ments that are given by 

/ _1 _1 I I \ 



m 



3 3 3 3 

0-10 



-10 3 2 
\ 1-10/ 

Quite generally, in the symplectic basis the periods have the expansion 



(3.32) 



n C0) = 2^ CfcPfcV' 



fc-i 



(3.33) 



fc=i 



where the cj, are complex constants and are constant four-vectors with components that are 
rational linear combinations of powers of a. 

As a further computational tools, we note in passing that for this model the Meijer periods 
(we find them following method developed in [2H] for quintic in the context of black hole 
attractors) are given by the following expression in Mathematica notations 



U 


(*) 


= MeijerG 




(*) 


= MeijerG 


u 2 


(*) 


= MeijerG 


U 3 


(*) 


= MeijerG 



{}j,{{0},{0,0,0}},-z 
{}},{{0,0},{0,0}},^ 
{}}, {{0,0,0}, {0}},-^ 
{}), {{0,0,0,0}, {}},z 



12 4 5 
6' 6' 6' 6 
12 4 5 
6' 6' 6' 6 
12 4 5 
6' 6' 6' 6 
12 4 5 
6' 6' 6' 6 

where U~ are valid for Im(z) < 0. 

The transformation to the Picard-Fuchs basis is given by the matrix 

/ 3 -5 3 -6 \ 

-14 3 

-10 

\ 3 -4 3-3/ 



so that 



V M^) J 



L ■ 



( Uo(^) \ 



with U~ 



for Im 



< 0. 



(3.34) 
(3.35) 
(3.36) 
(3.37) 



(3.38) 



(3.39) 
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3.3.1 Two vacua for one flux and an interpolating domain wall 

We now proceed to specifying the vacua for given values of the fluxes. Let us take, as a 
convenient choice among the solutions of the F-flatness conditions, the flux components to be 

h = {-50,-11,6,-4} (3.40) 

and 

f = -h-A 2 = {-4,-7,7,16}. (3.41) 

The virtue of this set of numbers is that the consequent vacua fall well within the conver- 
gence region of the Picard-Fuchs periods, thus making the numerical computations for the flow 
somewhat easier. 

It turns out that these fluxes give rise to two vacua: the first vacuum is a maximum at the 
Landau-Ginzburg point, with 

Vn = o, n = l±£fi . (3.42) 

The critical value of the potential is V(ipx,Ti) = —3e K \W\ 2 ~ —1219.5. We find the following 
diagonalized mass matrix for this vacuum 

mass = diag(-494.6, -494.6, -172.0, -172.0). (3.43) 
The second vacuum is a saddle point close to the conifold point, with 

tp 2 ^ 0.953 + 0.025i, r 2 ~ 0.162 + 1.088i, (3.44) 
and with potential equal to V(ip2, T 2) — —1269. In this case the diagonalized mass matrix is 

mass = diag(1260.5, -377.6, -341.3, 171.3). (3.45) 

The interpolating solution going from one vacuum to the other is found by solving the flow 
equations ()2.29j) . that become 

^ = -g^d^K + HW)) , (3.46) 
dr 

= -g^ df (K + MW)) , (3.47) 

du 

where we have used for convenience the variable u = In (a). 

The resulting flows are displayed below on the ip complex plane (Figure EJ) and on the r 
complex plane (Figure EJ). The corresponding flow of the potential V(if),r) as a function of 
u = In (a) is shown in Fig. [7| It shows clearly that there is no barrier. These domain walls 
are associated to flows that go from a maximum in the Landau-Ginzburg point, towards the 
conifold which is a saddle point. These are explicit examples of critical points which may be 
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Figure 5: Flow on ip complex plane. 
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Figure 6: Flow on r complex plane. 
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Figure 7: Flow of the potential V(ip(u), t(u)). 



sometimes repellers (for saddle points or maxima of the potential) and sometimes attractors, 
in particular stable attractors in case of minima of the potential. In this case, the flow tends 
to escape from the LG point towards the conifold. 
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We are confident that these results on interpolating domain walls are common to other 
one-parametric Calabi-Yau models that have been studied in |2Zj and generically have similar 
distribution of vacua on moduli space. 

4 BPS domain walls between vacua separated by barri- 
ers 

When the vacua are separated by a barrier, one can expect to find an absolute minimum as well 
as a meta-stable local minimum. This class of models in the stringy landscape was suggested 
in in addition to the AdS minimum and to the asymptotically Minkowski vacuum, there 
is also a supersymmetric Minkowski minimum at some finite value of the volume modulus. The 
main motivation for the new class of models with such a minimum was to have more flexibility 
with the gravitino mass, and to avoid certain cosmological problems specific to the standard 
KKLT scenario. Here we will show that the same class of models (which we will call KL models) 
allows two AdS minima, both at finite value of the volume modulus. Moreover, we will find two 
types of flows here: depending on the choice of parameters the superpotential may never vanish 
between the two AdS vacua, or it may vanish at some point and therefore change the sign along 
the flow. This results in several different types of BPS domain walls in the KL models. 

4.1 KL model with Minkowski and AdS minima 

The model has one modulus as in the KKLT case, but a more involved superpotential of the 
racetrack type |13| . The superpotential is a modification of the KKLT model by an additional 
exponential term 

W = W + Ae- ap + Be~ bp . (4.1) 

This model is a simple instance of multiple basins of attraction. It has one supersymmetric 
Minkowski vacuum at fixed volume, one supersymmetric AdS vacuum also at fixed volume 
and one supersymmetric Minkowski asymptotic vacuum at infinite volume as shown in an 
upper curve in Figure |H1 The plot of the potential clearly shows the barrier between the left 
Minkowski minimum and the AdS minimum, the superpotential is negative everywhere away 
from Minkowski vacua. For the choice A = 1, B = —2, Wo = —0.125, a = 7r/100 and b = tt/50, 

2 i 

also using once more the definition p = + ia, the vacua are at a = 0, 0a/i ~ 3.3 and 
4>AdS ~ 3.94, with V ~ —2.7 ■ 10~ 9 . Therefore, one can construct two interpolating domain 
walls. One is between the AdS minimum and the local Minkowski minimum, which we will call 
Mi. Another wall is between the AdS minimum and the Minkowski extremum at infinity, Moo. 
Here we study only the domain walls interpolating between Mi and AdS, see Figure 01 The 
features of the domain wall between M m and AdS minimum in this model will be the same as 
the ones studied in KKLT case in Section 13.11 In Figure El we plot on the top left the scalar 
potential V, and on the right the covariantly holomorphic superpotential Z, which vanishes at 
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Figure 8: The potential energy density V, blue upper curve, and the covariantly holomorphic super- 
potential Z, the red curve below, as the functions of the modulus (j) for the KL model. The potential 
is shown in units 10 -8 , and Z is in units 10~ 4 . Both V and Z have extrema at the same point, which 
is a reflection of the fact that supersymmetry is unbroken in the minima of V. This plot is for A = 1, 
B = -2, W = -0.125, a = tt/100 and b = vr/50. 
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Figure 9: KL wall between Mi and AdS. We plot the scalar <tJ>, Z (in units 10 -5 ), warp factor a and 
the curvature scalar R (in units 10 -8 ) as functions of the coordinate r. The curvature R is equal to 
zero in the Minkowski vacuum, has a peak at the domain wall, and approaches a constant negative 
value in the AdS vacuum. 



Mi and stays negative all the way towards AdS. The warp factor on the bottom left starts as 
a constant at Mi and interpolates into an exponent of r near the AdS. Finally, the curvature 
jumps from zero at Mi to some positive value and falls down to a negative constant at AdS. 
It is interesting that the curvature of the wall grows significantly between zero at Mi before 
becoming negative at AdS. 
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4.2 The KL model with two different AdS minima 
4.2.1 BPS walls with critical tension 

Now we are going to consider the models with two AdS minima. This can be done, e.g., by 
taking the KL model with a supersymmetric Minkowski minimum Mi, and change any of the 
parameters of this model. In particular, as one can see from Fig lTUl one can obtain an AdS 
minimum from the Minkowski minimum by a slight modification of a single parameter in the 
previous model: we changed A = 1 to A = 0.96. 

Using once more the definition p = e^» + ia, the AdS vacua are at a = 0, </>o ~ 3.3 and 

0i ~ 3.9, with V 0.034 • 10" 8 , and V l 0.29 • 10" 8 for a = 0. The potential V (upper 

curve) and the covariantly holomorphic superpotential Z (lower curve) for the KL model as 
function of the modulus are shown in Figure ITU1 for the choice A = 0.96, B = —2, Wq = —0.125, 
a = 7r/100 and b = 7r/50. Z is always negative. 



V 




Figure 10: The potential V (upper curve, in units 10 -8 ) and the covariantly holomorphic superpo- 
tential Z (lower curve, in units 10 ) for the KL model, as function of the modulus. Notice that we 
took A = 0.96 instead of A = 1, all other parameters being the same as for the potential in Fig. El 



We show various features of the KL wall between the two AdS vacua in Fig. ^2 We plot the 
scalar and the covariantly holomorphic superpotential, which is a negative constant near the 
first AdS vacuum and interpolates towards another negative value when it approaches another 
AdS vacuum. The warp factor a starts as a decreasing exponent and then changes into a faster 
decreasing exponent of r near the second AdS. The curvature R jumps from some constant 
negative value near the first AdS vacuum, becomes positive near the domain wall, and then 
falls down to a negative constant near the second AdS vacuum. 

The tension of this BPS wall is given by Eq. ()2.42|) . 

a = ^(IV^ - \V n . (4.2) 

It is critical in a sense explained in jjj: The tension of such BPS brane corresponds to a limit 
of the CDL bubble solution ^2] with infinite radius. When such BPS wall exist there is no 
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Figure 11: KL wall with a critical tension between AdS\ and AdS2- We plot the scalar and the 
covariantly holomorphic superpotential Z (in units 10 -5 ), which is always negative, warp factor a, and 
the curvature scalar R (in units 10 -8 ), as a function of the coordinate r. The curvature R has a peak 
at the domain wall, and approaches different constant negative values in two different AdS minima. 

tunneling since these supersymmetric walls saturate the CDL bound. We will discuss this issue 
later when we will study tunneling, see Section 

4.2.2 BPS walls with super-critical tension 

In the previous subsection we considered a KL model with two AdS minima and the super- 
potential which did not change sign on the way from one minimum to another. Now we will 
construct a KL model with two AdS minima and superpotential vanishing at some point be- 
tween them (Fig. 12). This can be done by a small deviation of one of the parameters: instead 
of taking A = 1 as in original KL model, or A < 1, as in the previous subsection with two AdS 
vacua, it is sufficient to take A > 1. 

For example, for the choice A = 1.05, B = —2, W = —0.125, a = n/100 and b = 7r/50, 

2 i 

also using once more the definition p = + ia, one finds that the AdS vacua are at a = 0, 

0o ~ 3.2, with V 0.8 • 10~ 9 , and at a = 0, 0i ~ 4, with V x 2.6 • 10~ 9 . The covariantly 

holomorphic superpotential Z is positive near the AdS minimum at <po . It crosses zero at 
~ 3.45 and becomes negative when approaching the second AdS minimum at 0i ~ 4. It 
approaches zero at infinite values of the volume modulus. 

The corresponding BPS solution is shown in Fig. We plot on the top left the scalar 
field 0, and on the right the covariantly holomorphic superpotential, which changes the sign on 
the solution. The warp factor on the bottom left starts as a growing exponent near the first 
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Figure 12: The potential V (upper curve, in units 10 -8 ) and the covariantly holomorphic superpo- 
tential Z (lower curve, in units 10 ) for the KL model as function of the modulus. Z is positive at 
4>q and negative at <j>\. Here we took A = 1.05 instead of A = 1, all other parameters being the same 
as for the potential in Fig. |HJ 

AdS minimum and changes into a decreasing exponent of r near the second minimum. Finally, 
the curvature jumps from some constant negative value near the first AdS vacuum, becomes 
positive near the domain wall, and then falls down to a negative constant near the second AdS 
vacuum. As far as we know, this solution provides the first explicit example of a smooth BPS 
domain walls in d — 4 with the superpotential vanishing between the critical points. 
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Figure 13: A BPS wall between two AdS domains in the version of the KL model where Z changes 
the sign on the domain wall solution. We plot the scalar <f>, the covariantly holomorphic superpotential 
Z (in units 10 -5 ), the warp factor a, and the curvature R (in units 10 -8 ). Notice the characteristic 
peak of the scale factor at the point where Z changes sign. 
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The tension of this BPS wall is given by Eq. (J23HJ). 

- = ^(N 1/2 + N V2 )- (4.3) 

It is super-critical in the sense explained in [J]. The tension of such BPS walls exceeds the CDL 
bound and therefore the tunneling is super-suppressed. 



5 Domain walls, near-BPS bubbles, and vacuum decay 

5.1 Decay of an uplifted vacuum to a collapsing universe 

In previous sections we studied domains with unbroken supersymmetry separated by static BPS 
domain walls. In our world supersymmetry is broken, but if supersymmetry breaking is not 
very large as compared to other parameters, one may hope that our domain wall solutions will 
provide an approximate description of slowly moving near-BPS domain walls. Just as in the 
situation with the near-extremal black holes, one may hope the physics of the near-BPS walls 
can be under a much better theoretical control than the physics of generic domain walls. It 
would be especially interesting to see whether the BPS domain walls may play any role in the 
cosmological context, dividing the universe into different regions separated by static domain 
walls. However, from the point of view of the Friedmann cosmology, the universe with a negative 
cosmological constant A < should collapse to a cosmological singularity within a very short 
time t ~ | yV | 1//2 . Therefore in the usual cosmological context one can hardly describe coexisting 
AdS domains separated by static domain walls. 

One of the most interesting applications of the results obtained in the last section is the 
theory of vacuum decay and bubble formation. It is well known that dS vacua obtained in the 
KKLT construction by uplifting of AdS vacua are unstable with respect to the decay to the 
Minkowski vacuum corresponding to infinitely large values of the volume modulus 0, i.e. to a 
supersymmetric 10D vacuum 

According to [12], the decay probability is given by 

P(0) = e ~ B = e- sw+So , (5.1) 

where S((f)) is the Euclidean action for the tunneling trajectory, and So = S(4>o) is the Euclidean 
action for the initial metastable configuration <fi = <po- 

The Euclidean action calculated for the metastable vacuum dS solution = O is given by 

So = -^f- < . (5.2) 

This action has a simple sign-reversal relation to the entropy of de Sitter space So: 

So = -S = + 2 -^f- . (5.3) 
Vo 
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Therefore the decay time of the metastable dS vacuum tdecay ~ P 1 {4>) can be represented in 
the following way: 

tdec ay = e s ^ = t r e s ^ . (5.4) 

Here t r ~ e s ° is the so-called recurrence time. It was shown in j3] that if the decay of the 
metastable dS vacua occurs due to the tunneling through the barrier with a positive scalar 
potential, then S(<p) is always negative, and therefore the decay always happens during the 
time shorter than the recurrence time t r . This result directly applied to the simplest KKLT 
model where the tunneling occurs through the positive barrier separating the metastable dS 
vacuum and the supersymmetric 10D vacuum. However, the situation with the tunneling to 
the minima with a negative vacuum density remained less clear; see Refs. |3U| l3*T| 1312] for the 
recent discussion. In what follows we will examine this issue using the results obtained in the 
previous section. 

Our consideration will be based on the observation of Ref. [7j that in all cases when Z does 
not vanish across the domain wall, the domain wall solutions can be represented as a limiting 
configuration describing the wall of the CDL bubbles of an infinitely large radius. For such 
bubbles, the tunneling action is also infinitely large, and the vacuum decay is impossible. In 
our context, this fact is related to the supersymmetry of the different vacua and of the 
interpolating BPS wall solution. 

However, in realistic situations supersymmetry is broken. For example, in the KKLT con- 
struction the metastable dS vacuum state appears after the uplifting of a supersymmetric AdS 
vacuum. Uplifting can be achieved, e.g., by adding a D3 brane contribution, which results in a 
supersymmetry breaking in the dS vacuum. Similarly, in the KL model one must uplift a 
Minkowski vacuum or an AdS vacuum to a dS vacuum, thus breaking the supersymmetry. As 
we will see, this may lead to a relatively rapid decay of the uplifted dS vacuum due to creation 
of the bubbles describing a collapsing universe with negative vacuum energy corresponding to 
the AdS vacuum. 

Note that in terms of a canonically normalized volume modulus 0, which we will use here- 
after, the D3 brane contribution looks as (Je~ 2v ^^^, where C is some constant. (The canoni- 
cally normalized field, 1/2 (d<p) 2 , is obtained by dividing the field <fi used in the previous section 
by the factor \/2.) This contribution rapidly decreases at large <p. Therefore this term may 
significantly uplift the AdS (or Minkowski) minimum at small values of <p, but it makes a much 
smaller effect on the AdS minima at large <fi. This suggests that the transitions should typically 
occur towards larger values of the volume of the compactified space. 

5.1.1 Decay of an uplifted Minkowski vacuum 

First of all, consider the KL potential with the supersymmetric Minkowski vacuum, and uplift it 
by a tiny amount to make the vacuum energy there equal to the present cosmological constant, 
Vo ~ 10~ 120 , in Planck energy density units. This should make this dS vacuum metastable, and 
lead to a vacuum decay due to creation of bubbles with negative vacuum energy. To simplify 
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Figure 14: A slightly uplifted KL potential, shown in units 10 
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the investigation of its decay, we will consider potentials such that the difference between the 
depth of the two minima is much smaller than the height of the barrier. An example of such 
a potential is provided by the KL model with parameters A = 0.1, B = —10, Wq = —0.00025, 
a = 27r/50 and b = 27r/25. For these values of parameters the potential has a supersymmetric 
Minkowski minimum at ~ 4.6 and an AdS minimum at <p ~ 5. If one slightly changes the 
parameters, e.g. takes A = 0.095, then the Minkowski minimum becomes an AdS minimum. If 
one uplifts the potential by adding a term Ce~ 2 ^^^, the potential acquires a dS minimum, as 
shown in Fig. El During all of these changes, the second AdS minimum changes only slightly. 

One expects that the thin wall approximation should be valid for the description of tunneling 
in such potentials. One part of this condition is automatically satisfied for the near-BPS 
bubbles. Indeed, by considering sufficiently small uplifting one can always make such bubbles 
arbitrarily large without changing the thickness of the wall. However, in order to use the thin- 
wall approximation as formulated in Ref . ^2] , one should also check that the difference between 
the depth of the two minima is much smaller than the height of the barrier, and that the width 
of the wall of the CDL bubble remains much smaller than the radius of the bubble even if 
one turns off the gravitational interaction. We checked numerically that this condition is also 
satisfied for our potential. 

Having in mind potentials of this type, we will use the thin-wall approximation to investigate 
the decay of the dS vacuum with the field O an d vacuum energy density V = V((f>o ) > to 
the state with the field (pi with a negative vacuum energy density V\ = V(cj) ) < 0. We will 
be interested in the realistic situation |V | <C |Vi|, as illustrated by Fig. H3I Examples of the 
vacuum decay considered in [T2] described transitions between dS and Minkowski space, and 
decay of the Minkowski vacuum. We will use the results obtained by Parke jH3], who generalized 
the results of ^2] for the tunneling between the states with arbitrary values of Vi in the thin-wall 
approximation. 

According to [54] . the bubble size is given by 



i i 9-^2_ i (Po_\ 
1 + ^ 4A 2 "I" { 2 a) 




4 ' 



(5.5) 
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where 

3a . 

Po = v^v, (5 " 6) 

is the bubble radius in the thin wall approximation ignoring gravity, a is the bubble wall tension, 
and 

A2 = v^v, (5 - 7) 

(it can be positive or negative; in our case it is negative), and 

A2 ^- < 5 - 8 > 

In the thin-wall approximation, the tension of the wall can be approximated by the integral 

a= / V(0) - Wo)] 1/2 • (5.9) 

J O 

Here 4>*i is the point near where V(0) — V(<po) = 0. (Ref. [T2] uses a slightly different 
prescription, which is equivalent to ours in the limit when the height of the barrier is much 
greater than the difference between V^.) In the limit when the radius of the CDL bubble becomes 
infinitely large, the wall tension is known exactly: 

- = ^(\Vi\ 1/2 ~ \Vo\ 1/2 ) , (5.10) 

see 0131], and also our derivation of this result in Sect. 2, Eq. (|2.42|) . 
Using the variables 

X= V2AJ ' V = W (5 - 11} 
one can represent Eq. (J5.5j) in the following convenient form: 



P 2 = o- (5-12) 

H l + 2xy + x 2 K J 

For a particular case of tunneling from Minkowski space (Vo = 0), the infinite bubble size regime 
(the threshold regime for the tunneling) corresponds to x — 1, y = — 1, when the denominator 
vanishes. 

The tunneling probability is given by e~ B , where 

B = 2(^$r <*■»>■ (5 - l3) 

Here the first term is the no-gravity result, and the gravitational correction factor r(x,y) is 
given by [HI] 

, 2[(1 + xy) - (1 + 2xy + x 2 )^} 
r(x,y) = — — f- 1 . (5.14) 

x 2 (y 2 - + 2xy + x 2 )^ 

32 



This factor also blows up at x — 1, y — — 1, i.e. the tunneling becomes forbidden. 

Now let us uplift the Minkowski vacuum, assuming that the AdS vacuum at <f>\ remains 
approximately at the same height V±. As we already argued, this is indeed the case if the AdS 
vacuum is at large values of the volume modulus 0. Uplifting changes not only Vo, but also 
the bubble wall tension a. Typically uplifting decreases the wall tension, which increases the 
probability of the tunneling; a sufficiently large uplifting may remove the wall altogether, see 

m 

Any small uplifting makes the tunneling from O to <f>\ possible. In order to study the size 
of the bubbles and the probability of the tunneling, let us change x and y by Ax and Ay. In 
the linear approximation Ax ~ Ay and both depend on the small uplifting Vo <C 1. In general, 
one would expect that the result should depend both on Ax and Ay. However, we find that 
in the leading approximation in Ax and Ay, the results of the expansion near x = 1, y = —1 
depend only on Ay, i.e. they are not sensitive to the small changes of the bubble wall tension. 
In particular, the denominator in (|5.12|) is equal to 2 Ay, where Ay = — for [ Vq | <C |Vi|. 
This yields 

Here Hi is the Hubble constant in dS space with vacuum energy density V . Thus, instead of 
an infinite size bubble, we have a bubble with the size equal to the size of the dS horizon! 

Now we will perform a similar investigation for the tunneling rate. By considering the 
function r~ 1 (x,y) in Eq. 1)5.14)1 and expanding it with respect to Ax and Ay, one finds that in 
the limit \V \ < |Vi| 

12tt 2 

B = —— = So/2 . (5.16) 

where So is the entropy of dS space, So = This means that the tunneling probability is 
given by the universal model-independent equation which has a simple geometric interpretation 
in terms of dS entropy: 

P~e- So/2 . (5.17) 

Several comments are in order. First of all, we have found that the decay rate of the vacuum 
obtained by uplifting of a supersymmetric Minkowski minimum does not depend on the tension 
of the wall, or on the depth of the supersymmetric AdS minimum. The only thing we needed 
to know is that prior to the uplifting these two minima were supersymmetric, which is a generic 
property of all known flux vacua. An unexpected feature of this result is that the tunneling is 
suppressed not by a factor e~ So , but by a factor e~ s °/ 2 . This means that the decay time of an 
uplifted Minkowski vacuum will be shorter than the recurrence time t r ~ e So by a huge factor 
e s °/ 2 ~ e io 120 /2 w m re turn to the discussion of this fact in Section |U1 
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5.1.2 Decay of an uplifted AdS vacuum 



As we have seen in the previous subsection, a tiny supersymmetry breaking achieved by uplifting 
of the supersymmetric Minkowski vacuum to the dS vacuum with Vo ~ 1CT 120 renders this 
vacuum unstable. Now we will extend our results for a more general case. 

Suppose that initially we have two supersymmetric AdS vacua, with depth Vb, V\ < 0, 
such that | Vo| < |Vi|. One can achieve it e.g. by considering the KL potential discussed in 
the previous section, but taking A = 0.09 instead of A = 0.095. Tunneling between the two 
supersymmetric AdS minima is forbidden. Let us also assume that the deeper minimum is at 
larger values of the volume modulus, so the KKLT uplifting of the minimum Vb to dS space with 
Vo ~ 10" 120 leaves the depth of the second AdS minimum almost unaffected. This is indeed the 
case in the model discussed above. We will also assume that prior to the uplifting |Vo| was much 
greater than the present value of the vacuum energy 10 -120 , so for all practical purposes it is 
not important whether we uplifted it to a dS state Vq ~ 10~ 120 or to the Minkowski state with 
V = 0. This subtle difference was important in our previous investigation because the degree 
of supersymmetry breaking was proportional to Vq ~ 10~ 120 . Now we are considering a much 
greater supersymmetry breaking, which is achieved due to a much more significant uplifting, 
from the initial value of Vb to (approximately) zero. In order to simplify notation and avoid 
mixing the values of Vo before and after the uplifting, we will assume that after the uplifting 
V = 0, and in all subsequent equations we will understand by Vo its value in the AdS minimum 
prior to the uplifting. 

In order to study the decay rate of the (nearly) Minkowski vacuum formed by uplifting, one 
should use, as before, eqs. ()5.13|) and (|5.14j) . In Minkowski vacuum, one has y = — 1, so the 
expression ()5.14|) looks singular, but in fact this singularity is fictitious, which can be easily 
seen by expanding it near y = — 1. 

The outcome of the calculations can be represented in the following form: 



Here a is the bubble wall tension for the tunneling in the uplifted potential. The first term 
in this expression is the probability of the decay of the uplifted Minkowski vacuum ignoring 
gravitational corrections. The second term represents the gravitational suppression of the 
probability of the tunneling due to the factor r(x,y). 

As we see, the probability of the tunneling depends on the wall tension a. However, the only 
part of the expression (j5.18|) where this dependence is crucially important is the gravitational 
suppression term (l — yjjfrp-) • If after uplifting we would still have two supersymmetric vacua 
(Minkowski and AdS), then we would have a = ; the last term in ()5.18|) would blow 

up, and the tunneling would be forbidden. However, prior to the uplifting we had a smaller 
value of tension, a = -^(IVi) 1 / 2 — (Vo) 1 / 2 ). As we argued before, uplifting tends to decrease the 
barrier between the two vacua, and therefore it tends to decrease the tension. (Our numerical 




(5.18) 
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investigation of the tunneling in the KL model with the potential shown in Fig. is consistent 
with this expectation.) If this is the case, the gravitational suppression gradually disappears, 
and the probability of the tunneling becomes very large. 

To illustrate various options, let us assume, in the first approximation, that a does not 
change during the uplifting. In this case we find that 



24tt 2 (C - 1) 
\Vo\ C 2 (2C-1) 



B = — ^ ' (5.19) 



where C 2 = > 1. 

For C = 0(1) (i.e. for V\ of the same order as V , prior to the uplifting), this expression 
coincides, up to a numerical coefficient g 2( 2 c-i) 2 = wr th the entropy of dS space with the 

vacuum energy | Vq | , see Eq. ()5.3|) . The main difference is that for the dS space in which we live 
now we have Vq ~ 1CT 120 , and the factor e _So is incredibly small, e _So ~ e -1 ° 12 °. Meanwhile the 
depth of the AdS vacuum prior to the uplifting can be a hundred orders of magnitude greater 
than 10~ 120 , and therefore the decay rate of this vacuum after the uplifting can be very large 
(see below for a numerical estimate). Note that if, as we expect, a slightly decreases during the 
uplifting in our model, then B will be even smaller, and the decay rate will be even higher. 

It would be nice to confirm our conclusions by other methods, or to calculate the tunneling 
probability directly, going beyond the thin-wall approximation. However, our main goal here 
was to show that generically there is no reason to expect a strong gravitational suppression of 
the decay of dS vacua if they are obtained by the KKLT uplifting of the supersymmetric AdS 
minima. Qualitatively, the main argument is as follows: If there are many AdS minima in the 
stringy landscape, and if the tunneling between them is suppressed by supersymmetry, then 
this suppression disappears after the uplifting and supersymmetry breaking. 

Let us check whether these intuitive expectations match our estimates, Eq. (|5.19|) . According 
to (|2.3jl . the scalar potential prior to the uplifting is given by V = \DiZ\ 2 — 3\Z\ 2 . In the 
supersymmetric minimum at 0o we have DiZ = 0, and Vq = —3\Z\ 2 . Uplifting practically does 
not change the position of the minimum <po, and therefore it practically does not change Z((po). 
Meanwhile the gravitino mass squared after the uplifting is given by m 2 ^ 2 = \Z\ 2 . This means 
that the depth of the AdS minimum before the uplifting is equal to — 3m 2 / 2 , where m 3 / 2 is 
the present value of the gravitino mass ^3] ■ Therefore our result for the tunneling rate of the 
uplifted AdS state (|5.19j) can be represented in the form directly related to the gravitino mass: 

8tt 2 (C-l) 4 

The scale of supersymmetry breaking in our world is much greater than the energy scale 
corresponding to the present value of the cosmological constant. For example, if the gravitino 
mass is m 3 / 2 > 10 2 GeV ~ 10~ 16 in Planck units, then m 2 ^ > 10~ 32 . Therefore Eq. (|5.2(J|) 
suggests that a typical decay rate of the AdS vacua after their uplifting to the dS state with 
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Vo ~ 10 is greater than e . Most importantly, this decay rate is greater than e ~ 
e -io 120 an enorm ously large factor e - 1034 / e - 10120 = e (io 120 -io 34 ) ^ e 10l2 °. 

This conclusion may have important cosmological implications which we are going to discuss 
in the next section. But before doing so, let us point out a possible caveat in our discussion. 
As we already mentioned, if the tunneling between the two vacua should go through the region 
where Z changes the sign, then a = -^(iVi) 1 / 2 + |Vo| 1//2 ). In this case, if we ignore the decrease 
of the tension during the uplifting, the tunneling between the two vacua does not occur. One 
can clearly see it by looking at the gravitational suppression term (l — jjr^-nr) • This term for 
a = ^(|Vi 1 1 I 2 + | Vol 1 / 2 ) is negative, which corresponds to super-suppression of the tunneling 

m 

On the other hand, we expect that in the realistic theories with hundreds of moduli fields 
there should be many AdS minima, some of which will remain AdS after the uplifting. Therefore 
it is natural to expect that the decay channels at least to some of these vacua will be open, 
and will be described by Eqs. (I5.19J) . ([5.20)1 . This is the only thing that matters to us, since 
the total rate of the vacuum decay is determined by the path of the least resistance, i.e. by the 
channels which maximize the tunneling probability. 



6 Sinks in the landscape and the wave function of the 
universe 

6.1 Tunneling, probabilities, and the wave function 

One of the main reasons to study vacuum transitions which occur on an incredibly large 
timescale is to find a probability distribution to live in one of the many vacuum states. This is- 
sue has a long history involving quantum cosmology, eternal inflation, and statistical properties 
of the stringy landscape. 

For a long time, vacuum transitions were studied in the context of inflationary cosmology. 
The main subject was the possibility of transitions between dS vacuum and Minkowski vacuum, 
or between two dS vacua. One of the most surprising and controversial results was obtained by 
Hawking and Moss [HS]. They found that a phase transition from a dS minimum with energy 
density Vq may occur via a quantum jump to the top of the barrier separating the two vacua, 
with the height V\ > V , and by a subsequent classical rolling to the second dS minimum with 
energy density V 3 . The probability of this jump, according to [351, was gi yen by 

p = e S 1+ S = exp + ( 6-1 ) 

\ vi J 

Initial interpretation of this result was rather obscure because the corresponding instanton 
seemed to describe a homogeneous tunneling, which was certainly impossible in an infinite (or 
exponentially large) inflationary universe. Moreover, from the derivation of this result it was 
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not clear why the tunneling should occur to the top of the potential instead of going directly 
to the second dS minimum. 



A proper interpretation of the Hawking-Moss tunneling was achieved only after the devel- 
opment of the stochastic approach to inflation j^Hl E3 EE EH]- 3 One may consider quantum 
fluctuations of a light scalar field with m 2 = V" <C H 2 = V/3. During each time interval 
St = H^ 1 this scalar field experiences quantum jumps with the wavelength ~ H^ 1 and with a 
typical amplitude S<p = H/2%. Then the wavelength of these fluctuations grows exponentially 
As a result, quantum fluctuations lead to a local change of the amplitude of the field 0, which 
looks homogeneous on the horizon scale H~ l . From the point of view of a local observer, this 
process looks like a Brownian motion of the homogeneous scalar field. If the potential has a 
dS minimum at 0o with m <C H, then eventually the probability distribution to find the field 
with the value at a given point becomes time- independent, 



That is why the probability to gradually climb to the local maximum of the potential at = 0i 
and then fall to another dS minimum is given by Eq. (jfi.lj) [3Q E3 EM EH] • 

Note that the distribution P(4>) is a fraction of the comoving volume of the universe cor- 
responding to each of the dS vacua. This probability distribution does not take into account 
different rate of growth of different part of the universe. Its interpretation can be understood as 
follows: At some initial moment one divides the universe into many domains of the same size, 
assigns one point to each domain, and follows the subsequent distribution P(<f>) of the points 
where the scalar field takes the value 0. 

A necessary condition for the derivation of this result in [3T"1 13*7] EBJ EHj was the requirement 
that m 2 = V" <C H 2 = V/3. This requirement is violated for all known scalar fields at the 
present (post-inflationary) stage of the evolution of the universe. However, one may try to 
generalize it for the situations with m 2 3> H 2 . Following [f""TJ (see also ""f""J HJ3 E] ) , we will look 
for a probability distribution Pj to find a given point in a state with the vacuum energy t^. The 
main idea is to consider CDL tunneling between two dS vacua, with vacuum energies Vq and 
V±, such that Vq < Vi, and to study the possibility to tunneling in both directions, from Vq to 
V\, or vice versa. 

The action on the tunneling trajectory, 5(0), does not depend on the direction in which 
the tunneling occurs, but the tunneling probability does depend on it. It is given by e~ s ^ +So 
on the way up, and by e~ s ^ +Sl on the way down |41j . Let us assume that the universe is in 
a stationary state, such that the comoving volume of the parts of the universe going upwards 
is balanced by the comoving volume of the parts going down. This can be expressed by the 
detailed balance equation 




(6.2) 



P oe -sw+s =Pie -sw +Sl 



(6.3) 



3 For a Hamiltonian approach to the HM tunneling see |4(J| . 
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which yields (compare with Eq. (jSHJ) 



Pi 0,0 / 24vr 2 24tt 2 \ 

e -5 1+ 5 = exp + ? (g_ 4 ) 



Po \ V V 1 

independently of the tunneling action S(4>). A 

Equations ()6.2|) and (|6.4|) imply that the fraction of the comoving volume of the universe in a 
state (ft (or 0i) different from the ground state 4> (which is the state with the lowest, but positive, 
vacuum energy density) is proportional to Co exp(|^y), with the normalization coefficient 

Co = exp(— ^fr-). The probability distribution Co exp ff^yl coincides with the square of the 
Hartle-Hawking wave function describing the ground state of the universe jH]. It has a simple 
physical meaning: The universe wants to be in the ground state 4>o with the lowest possible 
value of V(4>), and the probability of the deviations from the ground state is exponentially 
suppressed. This probability distribution also has a nice thermodynamic interpretation in 
terms of dS entropy S (IB] : 

P e Sl ~ s ° = e AS . (6.5) 



Po 

Here, as before, Si = — Si. This result and its thermodynamic interpretation played a substan- 
tial role in the discussion of the string theory landscape jB] . 

Unfortunately, there are some problems associated with this result. It suggests that the 
universe is similar to a pond filled with still water, where all things accumulate at the bot- 
tom: The fraction of the comoving volume with the smallest possible vacuum energy must 
be overwhelmingly large, P ~ exp ^^y-J- For a discussion of paradoxes associated with this 
conclusion and their possible resolution in the KKLT scenario see, e.g. [HH HH E] . 

Investigation of the stationary probability distribution alone does not give us a full picture. 
For example, the probability distribution (J6.2j) tells us about the fraction of the comoving 
volume of the universe in a given state, but it tells us nothing about the evolution towards 
this state. A partial answer to this question can be given by investigation of the stochastic 
diffusion equations describing the evolution of the scalar field in the inflationary universe. But 
now, instead of looking for the most probable outcome of the evolution, one should follow the 
evolution backwards and look for the initial condition <j) for the trajectories which bring the 
field to its final destination 0. In the stationary regime considered above, the corresponding 
solution looks very similar to (|6.2j) [39j: 

f 24vr 2 24tt 2 \ , . 

P{</>) -exp -—— + —— . (6.6) 



V{<t>o) V(<t>)J 

In this equation, however, O is n °t the position of the ground state, but a position of an arbi- 
trary initial point for the diffusion process which eventually brings us to the point <fi. As we see, 
the probability is maximized by the largest possible value of V(4>q). Interestingly, expression 

4 Some problems arising in the derivation of this result in the Coleman-De Luccia approach will be discussed 
in |H|. 
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exp(— y(^ ) ) describing the probability of initial conditions coincides with the expression for the 
square of the tunneling wave function describing creation of a closed dS universe "from noth- 
ing" |48j . whereas the second term looks like the square of the Hartle-Hawking wave function 
describing the ground state of the universe. In the stationary regime the squares of these two 
wave functions coexist in the same equation, but they provide answers to different questions. 

Here we should note that it is far from being clear whether the simplest probability distribu- 
tion which we studied so far should be used in anthropic considerations. Indeed, the probability 
that something happen at a given point, which is described by the probability distribution in 
the comoving volume as discussed above, does not take into account the exponential growth of 
the physical volume of the universe, which is the main property of inflationary models. For the 
discussion of the volume- weighted probability measure, see e.g. pM H9l lotH loT] 

Despite this reservation, it is quite important to study the probability distribution in the 
comoving volume since its investigation is an essential step in several different attempts to 
calculate the probabilities in the landscape. 

Here we would like to discuss another problem with the stationary probability distribution 
described above: This stationary distribution is not very useful during inflation when one can 
ignore the existence of the lowest dS space with Vo ~ 10~ 120 ; in order to obtain a stationary 
distribution one should take into account the growth of the physical volume of the universe 
39J. Moreover, as we are going to show now, this distribution does not necessarily apply to 
the string theory landscape, simply because in the KKLT scenario there are no stable dS vacua 
that could serve as a ground state of the universe. Metastability of dS space in the KKLT 
construction was emphasized in [3] and in many subsequent papers. Here we would like to look 
at this issue in a more detailed way. 

6.2 Probability currents and sinks in the landscape 

Stationarity of the probability distribution (j6.5J) was achieved because the lowest dS state did 
not have any further way to fall. Meanwhile, in string theory it is always possible for the 
metastable dS state to decay. It is important that if it decays by production of the bubbles of 
10D Minkowski space [5], or by production of bubbles containing a collapsing open universe with 
a negative cosmological constant |12j . the standard mechanism of return back to the original dS 
state does not operate any longer. 5 These processes work like sinks for the flow of probability 
in the landscape. The fraction of the comoving volume in dS vacua will decrease in time. 

To get a visual understanding of this process, assume that we are looking at the universe 
in the comoving coordinates, and paint black all of its parts corresponding to one of the two 
dS states, and paint white the parts in another dS state. Then the stationary regime which we 
discussed so far would imply that eventually the whole universe in average becomes gray, and 

5 One may speculate about the possibility of quantum jumps from Minkowski space to dS space [SSj, or even 
about the possibility of jumps back through the cosmological singularity inside each of the bubbles, but we will 
not discuss these options here. 
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the level of gray asymptotically will remain constant. 

Suppose now that some parts of the universe may tunnel to a state with a negative cosmo- 
logical constant. These parts will collapse, so they will not return to the initial dS vacua. If we 
paint such parts red, then the universe, instead of reaching the constant gray level, eventually 
will look completely red. 

This does not mean that eventually all parts of our universe will end up in the sink. Due to 
the exponential expansion of the universe in various dS vacua, the universe eternally rejuvenates, 
and the total volume of the universe in different dS states continues to grow exponentially 
PJ ESI IMl I55j . However, if one wants to find the ratio of the comoving volume of the universe 
in different states, then instead of solving the detailed balance equation (|6.3|) . one should solve 
the "vacuum dynamics" equations 

Pq = — Jo s — Joi + Jio , (6-7) 

Pi = -Jis - Jio + Joi ■ (6.8) 

Here Jq s = Poe~ Co is the probability current from the lower dS vacuum to the sink, i.e. to a 
collapsing universe, or to a Minkowski vacuum, J ls = Pi e 1 is the probability current from 
the upper dS vacuum to the sink, J i = Po e _So+ ' 5 ^'" is the probability current from the lower 
dS vacuum to the upper dS vacuum, and Jio = Pi e~ Sl+ ' 5( ^' is the probability current from the 
upper dS vacuum to the lower dS vacuum. Combining this all together, gives us the following 
set of equations for the probability distributions: 

P = _p e - Co - P e - s °+l + Pi e -Si+IW)l , ( 6 .9) 

A = _p ie -Ci _p ie -Si+IW)[ + p oe -So+|5(0)| (61Q) 

Since the entropy is inversely proportional to the energy density, the entropy of the lower level 
is higher, S > Si. Since the tunneling is exponentially suppressed, we have Si > so 
we have a hierarchy So > Si > \S(<p)\. We will associate the lower vacuum with our present 
vacuum state, with So ~ 10 120 , and therefore in what follows we will often have in mind the 
regime S > Si > \S(4>)\. 

For simplicity, we will study here the possibility that only the lower vacuum can tunnel to 
the sink, i.e. we will take the limit C\ — > oo and drop the term —J\ s = —Pi e~ Cl in Eq. (j6.10|) . 
We will analyze more general solutions of equations ()6.9|) . (j6.10|) in a separate publication |4*lj . 
On the other hand, we will keep in mind the results of the previous Section, where we have 
found that typically the probability of the decay of a metastable dS vacuum to a sink can be 
quite high, e~ c ° ~ exp(— 0(m~J 2 )) ^> e~ s ° ~ e _1 ° 120 . Therefore we expect that C -C S . 
Other mechanisms of a relatively fast vacuum decay have been discussed e.g. in j^Hl EH EE]- 
The main conclusions to be obtained in this section will be valid for any of these mechanisms. 
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By solving equations (jfi.lOJ) . one can show that P and -Pi exponentially decrease 

in the course of time due to the existence of the sink (the term — J 0s = —P e~ c °), but the 
corresponding physical volume of the universe exponentially grows, and the ratio Pi(t)/Po(t) 
approaches a stationary regime Pi(t)/Po(t) = p = const. In order to find p, one can add to 
each other our equations (without the term —Pi e~ Cl ). This yields 

(l+p)P = -P e- Co . (6.11) 

The solution is 

Pi / e- Co \ 

Po = — = -Po(t = 0) exp -— — t . (6.12) 

p V 1 +p ) 

As for the (asymptotically) constant ratio p = Pi(t) / P (t) , from eqs. (|6.9jh (j6.11|) one finds 
1 = (1+p) (1 + e co-s 0+ \sm _ pe Co-s 1+ \s(*)\) . ( 6 .i3) 

One may consider two interesting regimes. Suppose first that e~ Co <C e ~ Sl+ ' 5 ^-", i.e. the 
probability to fall to the sink from the lower vacuum is smaller than the probability of the 
decay of the upper vacuum. In this case one recovers the previous result, Eq. (|6.5)l . which is 
related to the square of the Hartle-Hawking wave function: 

p = § = e Sl ~ So < 1 . (6.14) 
Po 

Now let us consider an opposite regime, and assume that the decay rate of the uplifted 
dS vacuum to the sink is relatively large, e~ c ° ^> e _Sl+ ' s ^' (which automatically means that 
e~ c ° > e ~ So+|swl ). In this case the solution of Eq. flFTTT^ is 

p = ^ = e s ^ s ^~ Co w e Sl "l s WI » 1 , (6.15) 

i.e. one has an inverted probability distribution. This result has a simple interpretation: if the 
"thermal exchange" between the two dS vacua occurs very slowly as compared to the rate of 
the decay of the lower dS vacuum, then the main fraction of the volume of the dS vacua will 
be in the state with higher energy density. 

Note that the total fraction of the comoving volume remaining in both of the dS spaces 
exponentially decreases in time due to the existence of the sink in the lower dS vacuum: 

Po + P^expf-^-^t) . (6.16) 

This result may have interesting methodological implications. In many recent discussions of the 
probabilities in the landscape the authors consider the world as seen by an imaginary eternal 
observer. In the context of the stringy landscape such considerations are problematic, because 
every would-be eternal observer living in the landscape (unlike an observer in an eternally 
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existing dS space) begins his life in a cosmological singularity and ends up in a sink, within the 
typical time determined by Eq. (jfi.lfij) . Moreover, if the sinks occur in the higher dS vacua (see 
below), then the "eternal observers", or markers, will have a good chance to die well before 
reaching the thermalized regions produced after inflation. Introduction of a finite lifetime of 
such observers may result in the modification of the probability measure as suggested in [5T] . 

Here we should make an important clarification. The probability distribution Pi describes 
the fraction of the volume of the universe in a particular dS state. However, when the bubbles 
of a new phase expand, their interior becomes an empty dS space. If we are usual observers 
which are born after reheating of the universe, then one may argue that the probability to live 
in the bubble dSi is proportional not to the relative volume of the bubbles p = Pi/Pq but to the 
frequency of their production, i.e. to the ratio of the probability currents Joi/^io jUIlED!- In the 
absence of the sink, the fraction of the comoving volume which flows to the lower dS vacuum 
due to the tunneling from the upper dS vacuum is equal to the fraction of the volume jumping 
upwards from the lowest vacuum to the higher vacuum. In other words, the two probability 
currents are exactly equal to each other, 

Joi = Jio , (6.17) 

which is the essence of the detailed balance equation (|6.3|) . Interestingly, our new results 
imply that this regime remains approximately valid even in the presence of the sink, under the 
condition e~ c ° < e - s i+\ s W\. 

On the other hand, in the regime described by Eq. (j6.15|) . which occurs if the decay rate to 
the sink is large enough, e~ Co ^> e~ Sl+ ' 5 ^', one has a completely different result: 

T P c -So+\S(4>)\ 

Jio p ie -si+|swi v i 

Thus we have a crucial regime change at the moment when the decay rate of the lower vacuum 
to the sink starts competing with the decay rate of the upper dS vacuum. 

Now we will compare the probability current to the sink, J 0s , and the probability current 
from the lower dS upwards, J i- in the absence of the sink, J 0s was zero, so every point falling 
to the lower dS vacuum eventually jumps up to the higher dS vacua. This leads to eternal 
recycling of dS vacua |32J |HJ 113 , and to the problems discussed in |13J H?j : In the absence of 
the sinks (i.e. in the models with eternal dS space) each point in the comoving volume would 
indefinitely wonder between different dS spaces. In this scenario, it would be more probable 
that our part of the universe was formed as a result of a large quantum fluctuation upwards 
from an empty dS space, instead of being formed at the stage of inflation. This would leave the 
homogeneity of our universe, and the origin of perturbations of metric unexplained [4*31 147] . 

Invention of the KKLT mechanism demonstrated that all dS vacua in string theory are 
metastable with the lifetime smaller than the recurrence time t r ~ e s ° ~ e 10l2 °. This means 
that we do not have much time for the incredibly improbable jumps upwards to happen |47j . 
However, the time necessary for the improbable fluctuations upwards, which could re-create 
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conditions for the existence of life, is also somewhat smaller than the recurrence time |47j . 
Therefore in order to make sure that the universe does not recycle back from our dS state we 
would need to double-check that the probability to jump upwards is indeed strongly suppressed 
as compared to the probability to tunnel to the sink. 

The estimates of the decay probability of dS vacuum contained in our paper confirm that 
this is indeed the case. Indeed, investigation of the tunneling in several models discussed in our 
paper suggests that the probability that any given part of dS space with V ~ 1CT 120 will jump 
upwards and participate in the process of recycling of the universe is suppressed by a factor of 
0(e~ 1Ql ). This can be expressed in terms of the ratio of the probability current from the lower 
dS vacuum to the upper vacuum Joi with the probability current from the lower dS vacuum to 
the sink J 0s : 

= e -So+\sw\+c w e -s _ e -io 120 _ / g lg x 

Jos 

The last part of this equation follows from the condition So ^ Si > 15(0)1, and from the 
estimate of the decay rate obtained in the previous Section, which suggests that in many cases 
this rate is quite high, Co So- 

Despite this fact, recycling may still be possible, since for e~ Co <C e~ Sl+ ' s ^' one has Joi = 
J\q. On the other hand, for e~ Co ^> e~ Sl+ ^ s ^ one has Joi — e _1 ° 12 ° J 10 . This suggests that in 
this regime the probability that our universe was created due to the jumps up [I3J EZ] becomes 
exponentially small as compared with the probability that it was formed as a result of inflation. 
Using our results, one can easily check that this is indeed the case if the rate of spontaneous 
formation of the parts of the universe of our type in the vacuum dSo is smaller than the decay 
rate of the vacuum dSo to the sink. For wide sinks, this condition can be easily satisfied. 

This conclusion is closely related to the inversion of the probability distribution for e~ Co ^> 
e -Si+|S(0)| (widg s i n k; regime), see Eq. (|bM5|) . Note that inverted probability distributions in 
the presence of the sink represented by Minkowski vacuum naturally appear in the simplest 
chaotic inflation scenario with the volume- weighted probability measures proposed in |89| I49j : 
see Figs. 10 and 12 in [39 , where the basic theory of the volume- weighted distributions in the 
inflationary cosmology was developed. With the first of the two volume-weighted probability 
measures proposed in (SHI EH] the problem of non- inflationary jumps up does not appear at all; 
the growth of the total volume of the universe (and the corresponding fluxes) is predominantly 
determined by the exponential expansion of domains with larger vacuum energy and by the 
subsequent inflation on the way down. 

In our investigation so far we considered the simplest model with two dS vacua, with the 
probability leak in the lower dS vacuum. On the other hand, if the probability leak occurs only 
in the upper dS vacuum, one may ignore the term — J 0s = —P e~ Co in Eq. (j(19|) but keep the 
term — Ji s = —P\e~ Cl in Eq. (jfi.lOj) . In this case, for e~ Cl <C e~ Sl+ ' s ^! one again recovers the 
distribution (jdl4j) . whereas for e~ Cl ^> e~ Sl+ \ s ^ one finds 

p = ^ = e-^+l^)!^ w e- So ~ e- 10l2 ° . (6.20) 
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Note that in this regime we have a flat probability distribution, which does not depend on 

More generally, one should consider a dynamical equilibrium of a system of many dS, 
Minkowski and AdS vacua, where the probability leaks may occur at each level. The frac- 
tion of the comoving volume occupied by different metastable vacua in such a universe will be 
exponentially sensitive to the decay modes for each of these vacua. This makes the investiga- 
tion of the probabilities in the landscape much more involved, but also much more interesting. 
Previously we were trying to construct the wave function of the ground state of the universe. 
Now we should learn how to study the universe with the holes in the ground. 

7 Discussion 

There are two possible approaches to string theory landscape. The first one is based on the 
calculation of the wave function of the universe, which may describe all branches of the wave 
function corresponding to all possible vacua. In fact, the first two proposals for the anthropic 
solution of the cosmological problem in the context of inflationary cosmology were based on the 
investigation of quantum creation of the universe with various values of scalar fields and fluxes 
[2H], or with different types of compactification [oTi] . 

Another approach is to consider an inflationary universe consisting of many parts corre- 
sponding to all possible vacua [T], which allows, in particular, to consider a single universe 
divided into parts corresponding to all possible values of scalar fields, fluxes and compactifica- 
tions, and look for the anthropic solution of the cosmological constant problem in this context, 
see e.g. jHU E21 E31 EH El EH] ■ Conceptually, this approach is much simpler than the approach 
based on quantum cosmology, but in order to make it sufficiently general one must find a way 
to describe the universe consisting of many different parts and describe possible transitions 
between them. This was the goal of our paper. 

One of the problems here is that a fully consistent "landscape action" with the bulk and 
brane sources is not known [HHl EH EH]- However, one can hope that such action can be 
constructed at least in the case of unbroken supersymmetry following the "supersymmetry in 
singular spaces" proposal [HH] • This strategy was already applied in examples of fully consistent 
super symmetric bulk and brane actions in [7U| which include the D8 domain walls with the 
piece-wise constant Go flux. In a more general situation with moduli stabilization when various 
fluxes may jump due to brane sources, the corresponding construction has still to be developed. 

Therefore in the current project we restricted ourselves to studies of the regions of the 
landscape where stabilization of the moduli is under control and fluxes are fixed. Even in this 
restricted and simplified situation we were able to find multiple vacua and BPS domain walls 
separating them. 

In the first part of this paper we concentrated on the stringy landscape prior to the uplifting 
and found many BPS domain wall solutions separating different parts of the stringy landscape 
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from each other. We did it using the methods developed earlier in [7], and also more recent 
methods based on the "new attractor equations" |17j . 

In the second part of the paper we studied what happens with the domain walls after the 
uplifting. In general, potentials of many moduli before and after uplifting may have many 
different AdS minima of the moduli potential, for the same values of fluxes. Before the uplift- 
ing, the tunneling between these minima was forbidden by supersymmetry. One could expect, 
therefore, that after the uplifting the rate of tunneling may be related to the degree of super- 
symmetry breaking. Indeed, in the particular examples that we have studied the probability 
of the decay of the uplifted vacuum to the remaining AdS vacuum was suppressed by a factor 
~ exp(— 0(m^ 2 )) , where 7713/2 is the gravitino mass. This decay probability is much greater 

j , _a _in 120 

than e ~ e . 

One of the interesting features of the decay of dS vacua to Minkowski or AdS vacua is that 
this decay is irreversible. This process is substantially different from the stationary "thermal 
exchange" between various dS vacua. The existence of the channels of the irreversible vacuum 
decay, which serve as sinks for the flow of probability, is a distinguishing new feature of the 
stringy landscape which deserves further investigation. As we show in Section 6 of our paper, 
the existence of the sinks strongly affects the probability distributions in string cosmology. 
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